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Summary. — The Gamma Ray Inverse Problem is discussed. Four methods of 
spectral deconvolution are presented here and applied to the BATSE's MER data 
type. We compare these to the Band spectra. 
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1. Introduction 

The Gamma Ray Inverse Problem ([!]) is as follows: from a given number of observed 
parameters we have to reconstruct the intrinsic spectra. This will be estimated at some 
other number of points. The relation between the measured parameters and the intrinsic 
spectra is given by the detector response matrix (DRM) . In the case of the BATSE MER 
data, the spectra has to be approximated at 62 intervals from 16 measured count rates. 



The inverse problem corresponds to the discretized version of the Fredholm integral 
equation of the first kind. 



(1) 



d{y) = I R[x,y)f{x)dx. 



Here R{x,y) is the so-called kernel function, in our case the detector's response matrix. 
d{y) is the measured count rate, f{x) is the quantity to be determined. This is known 
to be an ill-posed and underdetermined problem. In the case of BATSE MER d is a 16 
element vector, / a 62 element vector and i? is a 16 x 62 matrix. 
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2. - Data 

Data were taken from the BATSE database (ftp:/ /cossc.gsfc. nasa.gov). We only used 
the Medium Energy Resolution data type. Baekground fits were made using intervals 
by Norris (http:// cossc.gsfc.nasa.gov/ docs/ cgro/ batse/ batseburst/ sixtyfour_ms/ 
bat_files.revamp_join) for the cat64 data. We have subtracted 3'"'' degree polynomials 
from all the 16 channels of MER data. 

3. - Methods 

Throughout the literature there are articles which discuss in detail the inverse problem 
in astronomy ([1, 2]) and other sciences ([3]). Gamma-ray bursts' spectra are most 
commonly reconstructed by the method parameter fitting also knows as forward folding 
[4]. This is done by postulating spectral forms with some adjustable parameters. Once 
multiplied with the DRM we get a solution for the count rates and we then compare 
them to the measured set of data. We accept the set of parameters which minimizes the 

of the solution and the measured counts. Although widely used and a well established 
method, still it is obliging when it comes to the shape of the spectra and it consists of 
estimating the parameters and not the spectra itself [1]. 

3'1. Singular Value Decomposition . - We decompose tlie response matrix using this 
well known mathematical method. We solve the GRIP by miiltiplying with the gener- 
alised inverse obtained through SVD. Spectra obtained from simply multiplying with the 
general inverse tend to be noisy. To avoid this we can use filters. 

3'2. Maximum entropy method. - This method introduces some a priori information 
to reconstruct the spectra (m). This is done by fitting a Band GRB spectral form [5]. We 
propose to maximize the expression: L{f, A) = XS{f, m) — where A is a parameter, 
m is the a priori spectra (from model fitting or SVD). It is with respect to this that 
we measure the entropy of the intrinsic spectra (/). The expression of the entropy is 

as follows: S{f,m) = J^fli^"^^ fj ~ '^j ~ fj^^S, (m^)' When L has a maximum then 
VL(/, A) = 0.. Prom this we derive: Alog (^^) = R^[(T-'^]{d - Rf), where [a^] is a 
diagonal matrix with elements of the tr^ vector and R is the response matrix. Since 
log {j^ results from the linear combination of the i?^'s cohuxms, we can search for the 

solution in the form of: / = me^~^^. Here t« is a 16 element vector. This ensures that 
/ will always be positive, and since w has a smaller scale, it is easier to find it as a 
result of a minimizing process. The choice of A can be made with the so-called cross- 
validation technique or another method based on a Bayesian approach. We have applied 
the method to BATSE MER data. First results suggest that the Band GRB function 
provides a good enough fit. 

3'3. Backus Gilbert method. - Its main advantage is that it is a linear method([6]), 

data in the intrinsic spectra is reconstructed as a linear combination of the measured 
count rates: f{Ej) = YhLi ai{Ej)di = Y.k=i ^k{Ej)fk and Ak{Ej) = Y^iii ai{Ej)R^^k- 
We call R-iE) the resolution function. Our aim is to approximate the resolution function 
with Dirac-delta ftmction as well as possible. By doing this we minimize for a,j. The 
other issue at hand is the variance of the spectra. This is dealt with by minimizing 
<^j{Ej) = Ei=iai(^^j)f^i also for o^. These two expressions are combined into a single 
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Fig. 1. - Sample Backus Gilbert inverse. 

one depending on the parameter 9 in the following manner: w{Ej,9) = TZ{Ej) cos0 + 
va]{Ej) sine. By minimizing w with respect to we get the Backus Gilbert inverse 
(Fig.l). {v is merely a factor to get the two components to the same order of magnitude) 

3'4. Phillips- Towmey method. - We will provide only a brief description of the method. 
Like at the maximum entropy method, we are also in need here of an a priori spectrum. 
This information is used to reduce the characteristic length scale of the problem. We 
introduce the sum-of-squares of the spectrum's k-th derivative (Cfe). 

We propose to minimize some linear combination of the and the Ck- As an a priori 
spectrum we fit a Band GRB function. We find that there is only insignificant gain in 
the function to be minimized with respect to the Band GRB function. 

4. — Conclusion 

Having applied several of the above presented methods to invert BATSE MER spectra, 
wc conclude that the methods give a consistent estimation of the incident spectrum each 
with regard to the minimizing function they use. The SVD method is the fastest, the 
Backus Gilbert requires the most time, but it is less sensitive to noise. We have also 
found that the Band GRB function provides a good fit to the spectra we analyzed. 
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